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D. Peiia Arteaga/ M. Grasso/ E. Khan,^ and P. Ring^ 

^ Institut de Physique Nucleaire, Universite Paris-Sud, IN2P3-CNRS, F-91406 Or say Cedex, France 
"^ Physikdepartment, Technische Universitdt Miinchen, D-85748, Garching, Germany 

Covariant density functional theory is used to study the efTect of strong magnetic fields, up to the 
limit predicted for neutron stars (for magnetars B ~ 10^*G), on nuclear structure. All new terms 
in the equation of motion resulting from time reversal symmetry breaking by the magnetic field and 
the induced currents, as well as axial deformation, are taken into account in a self-consistent fashion. 
For nuclei in the iron region of the nuclear chart it is found that fields in the order of magnitude of 
lO^^G significantly affect bulk properties like masses and radii. 
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I. INTRODUCTION 

Several studies (e.g. [ll, [2']) have determined the pres- 
ence of intense magnetic fields, up to 10^^ G, on the sur- 
face of neutron stars. Theoretical models suggest that 
these magnetic fields might reach up to B w 10^^ G, and 
even larger values if one considers the limit imposed by 
the virial theorem (B w 2 • 10^® G 01). The infiuence of 
the magnetic fields on the equation of state (EOS) has 
been thoroughly studied and reported in last years (e.g. 
[S-Q). It should be noted lower than 10^® G magnetic 
fields can influence the low density parts of a neutron 
star, such as its surface layer. 

The outer crust is fundamentally composed of well sep- 
arated nuclides, and its structure determined by the en- 
ergies of isolated nuclei, the kinetic energy of electrons 
and the lattice energy. Thus its composition depends 
very much on the binding energy of stable and unstable 
nuclei in the outer crust below the neutron drip density 
and of neutron-rich nuclear systems above the neutron 
drip. At the lowest densities it is thought that the most 
abundant component is ^^Fe because of its high binding 
energy, a fact actually observed in emission spectra (e.g. 

0)- ^ 

The magnetic field strength required to directly influ- 
ence the EOS can be estimated by considering its effects 
on charged particles. Charge-neutral, beta-equilibrated 
neutron star matter contains both negatively charged 
leptons (electrons and muons) and protons. Magnetic 
flelds quantize the orbital motion (Landau quantization) 
of these charged particles. When the Fermi energy of 
the proton becomes signiflcantly affected by the mag- 
netic field, the composition of matter in beta equilibrium 
is modified. This is refiected in a change of the pressure 
of matter. It has been found in Ref. ^ that this occurs 
for fields of approximately 10^^ G, and that in general 
leads to a stiffening of the EOS. 

However, there are very few studies [3,[3] of the changes 
that these very intense magnetic fields may eventually 
bring to the composition of the crust. The structure and 
composition of the crust is important in the thermal and 
rotational evolution of neutron stars, in particular in the 
theory behind glitches [l|. Some other studies [ifl, [III 



point to a magnetically driven crust activity as the source 
of soft gamma repeaters (SGR). 

So far, the impact of intense magnetic fields in nuclei 
found in the outer crust has been studied on a qualitative 
way using a simple non self-consistent method by Kon- 
dratyev and collaborators [1, Q • It was found that fields 
as low as 10^^ G may modify the nuclear shell structure, 
well within the range of theoretically possible magnetic 
strengths. There are, however, still some questions that 
need to be addressed: i) What is the minimum field that 
is able to significantly alter the nuclear structure? ii) Is 
this field low enough to be found in a significant propor- 
tion of neutron stars or magnetars? iii) Is this effect big 
enough to influence astrophysically relevant situations 
and processes, e.g. neutron star outer crust composition 
or final element abundances in nucleosynthesis scenarios? 

The objective of the present work is to try to find an- 
swers to these questions, in a quantitative way if possi- 
ble, using a fully microscopical description of the nuclear 
system within covariant density functional theory. The 
formalism used will be introduced in Section II. A general 
discussion of the effects of the external magnetic field on 
nuclei will be given on Section III, and particularized to 
an example nucleus. A discussion of the possible infiu- 
ence of the magnetic fields on neutron star outer crust 
nuclei can be found on Section IV. Finally, Section V is 
devoted to the conclusions. 



II. FORMALISM 

Covariant density functional theory starts from an ef- 
fective Lagrangian, that includes the nucleon and as 
many meson fields as needed to reproduce basic nuclear 
properties like saturation (a detailed discussion can be 
found inRefs. [l3,[ll ): 

C = Cn + Cm + Cint + CbO + CbM- (1) 

Cn refers to the Lagrangian of the free nucleon 

Cn = i'itl^'d^ - m)^, (2) 

where m is the bare nucleon mass and ijj denotes the 
Dirac spinor. Cm is the Lagrangian of the free meson 



fields and the electromagnetic field generated by the pro- 
tons 

- \r^,.W- + \vfilp^(P - \f^,,F^^ + C/(a), (3) 

where mo-, r^lu)^ rUp are the meson masses, and U{a) = 
(52/3)17^ + (33/4)0-'* is the standard form for the non- 
linear coupling of the a meson field. The interaction 
Lagrangian Cmt is given by minimal coupling terms 



gpi)^^Tp^'^ - e'lp'^'^A^ip, 



(4) 



where g^, g^, gp and e are the respective coupling con- 
stants for the CT, cj, p and photon fields and, of course, 
e vanishes for neutrons. In the previous and subse- 
quent formulae, bold symbols denote vectors in ordinary 
space, and arrows vectors in isospin space. These three 
terms, Cn, J~-m and Cint compose the standard RMF 
Lagrangian. Throughout this work the parameter set 
NL3 [14l| is used for the masses and coupling constants 
of the model. This parameter set has been thoroughly 
used and has led to very successful description of many 
nuclear properties. In addition, there are two new terms 
corresponding to the interaction of the nuclear system 
with an external magnetic field: i) the coupling of the 
proton orbital motion with the external magnetic field. 



Cbo = -e^7^4^)V, 



(5) 



and ii) the coupling of both proton and neutron intrin- 
sic dipole magnetic moments with the external magnetic 
field il5|J 



where 



Here 



^ fii- 



CbM = -V-Xr^ V", 



xt:) - ^^r.mi^cj^^F^'^^^^r 



= §[7m'7i'] and p.N = eh/2r 



(6) 



(7) 



the 



clear magneton and k„ = 5„/2, Kp = gp/2 — 1 with 
gn = —3.8263 and gp = 5.5856 are the intrinsic magnetic 
moments of protons and neutrons. Interactions with the 
external magnetic field are marked the superscript (e). 
This field is considered to be externally generated, and 
therefore there is no associated field equation and thus 
no other bosonic terms in the Lagrangian. Both terms 
J~-BO (0) and Cbm & have to be taken into account since 
at the magnetic field strengths of interest {B w lO^^G) 
they are of the same order of magnitude. 

The Hamiltonian density can be derived from the La- 
grangian density of Eq. ([T|) as the (0,0) component of 
the energy-momentum tensor, leading the to the energy 



functional Eb[p,4'] (see in Ref. flS'l for details). 

Eb[p, 0] = Tr [{a (-^V - eA^^^) + /?(m + x^^J)) p 
-f^Tr[(/3r„0„)p] 

m 



(l9^0m)^ + TrTOm^m 



(8) 



where the upper sign holds for scalar and the lower sign 
for vector mesons. 



/5W = Ei'/''M)<v^^M 



(9) 



is the relativistic single-particle density matrix and the 
traces run over the Dirac indices and over the integral in 
r-pace and, according to the no-sea approximation, the 
index i runs over all the occupied levels in the Fermi sea. 
The index m — a,uj, p,e runs over the various meson and 
electromagnetic fields and the vertices F™ read 



.9<T, r^ - gui", 



Fi'-fl„f7^, r^: = e-f^, 



(10) 



and \ml^(t>l^ 



has to be replaced by im^cr^ -I- U{a) in the 
case of the tr-meson. It is customary at this point to in- 
troduce an additional term into the energy functional to 
account for pairing correlations [la|, at least in its sim- 
plest BCS approximation. In the present study, however, 
pairing effects shall be completely neglected. It is a well 
known fact [13] that static magnetic fields lead to a re- 
duction in pair correlations in superconductors, and to 
the appearance of a critical field where all such correla- 
tions vanish. 

The functional (|8]) follows the spirit of magnetic-field- 
and-density functional theory (BDFT) [l^, in which the 
vector potential is introduced as an explicit dependence 
in the energy functional. Considering that astrophysical 
magnetic fields can be taken as constant on a nuclear 
scale (i.e. their functional form is fixed), it would be 
of little advantage to use the more general current-and- 
density functional theory (CDFT) [1^ , which generalizes 
density functional theory (DFT) with the inclusion of 
an external vector potential in a universal fashion. And 
since there is no practical value in considering the ex- 
ternal magnetic field B as an independent variable, it 
shall be regarded in the density functional as a paramet- 
ric variable. Minimization with respect to the density p 
in the Hartree approximation |13| and considering only 
static configurations leads to the stationary Dirac equa- 
tion for the nucleons and to the Klein-Gordon equations 
for the mesons 

hoipi = Citpi, (11) 

[-A + mf„] (b,n^Tj2 ^^^rn1p^, (12) 



where m^„ has to be replaced by m^ + 520" + ffso"^ in the 
case of the cr-nieson and where the Dirac Hamihonian 
has the form 



5p 



= a (-zV - V) + Fo + /3(m + 5) + /3xi:\ (13) 

with the scalar and vector potentials S and V^ defined 
as 

S = -g^a, (14) 

Vfj, = g^ujf, + gpT3Pf,^3 + eAf, + eA^^'' , (15) 

At this point it may be useful to fix the functional form of 
the magnetic field. Of course, it is not constant through- 
out the neutron star. However, the scale of changes 
is much larger than the microscopic nuclear scale [3|. 
Thus, the magnetic field B within each individual nu- 
cleus might be considered constant. If one chooses the 
intrinsic z-axis in the direction of this constant external 
magnetic field B = (0,0,5) and cylindrical coordinates 
(z, r, if) the contribution of this external field to the vec- 
tor potential V can be written, in the symmetric gauge, 
as 



Ai-^) 



rB 



(16) 



where e^p is the unit vector associated with the azimuth 
angle (/3 and r is the distance from the symmetry axis. As 
discussed in Eq. (|4ip of the Appendix B the contribution 
of the intrinsic magnetic moments is given by 



At3 



-KraMJvSsS 



(17) 



It is clear that the presence of the magnetic field 
breaks spherical symmetry for the Dirac and Klein- 
Gordon equations. Only axial symmetry is preserved for 
fields of the form ([T5|) . As discussed in Ref . [l2] the spinor 
solutions in Eq. (jlip can be written, in axial symmetry, 
as 



\Mr)) = 



1 

2tt 



I /+(r,z)e^("'-i/2)^ \ 

/r(r,z)e*(f^-+i/2)^ 

i5+(r,z)e'(^^'-i/2)'^ 



V ^9^ ( 



r, z)e 



i(ni + l/2)i^ 



Xt.(t). (18) 



They are characterized by the angular momentum pro- 
jection 51, the parity tt and the isospin projection t. For 
even-even nuclei and in the absence of an external mag- 
netic field, according to Kramers rule, for each solution 
■0^ with positive 57^ there exists a time-reversed one with 
the same energy, denoted by a bar, 1 := {ei, —11^, 7r,i}. 
However, time-reversal symmetry is broken by the mag- 
netic field, so the two-fold degeneracy is not present and 
one needs to consider both solutions separately. This 
breaking of time-reversal symmetry in the intrinsic frame 
leads to the appearance of time-odd mean fields and non- 
vanishing currents which induce space-like components 



of the vector mesons w and p, usually referred as nuclear 
magnetism |20l - [23 |. It is a great advantage in relativis- 
tic nuclear density functionals that these time-odd mean 
fields are determined by the same coupling constants gi^ 
and gp as the well determined time-even fields. 

In non-relativistic nuclear density functionals such as 
Skyrme [2j, [2J| or Gogny y^ there are, in principle, 
also relations connecting time-even and time-odd parts 
through Galilean and gauge invariance [2J]. However, 
these relations do not connect spin and spatial degrees 
of freedom as consistently as Lorentz invariance and, in 
addition, there is ambiguity, because many of these very 
successful functionals, still in use, are adjusted without 
taking them into account. 

As shown in Ref. 7m\ for fields of the form (IT5|) all 
induced currents and magnetic potentials are parallel to 
e.^ and axial symmetry is preserved as a self-consistent 
symmetry [26]. One can write out explicitly the Klein- 
Gordon equations for the time- and space-like meson 
fields as 

(-A + m^) cT = -.9, {pl + p:) - g^o-" - g^a\ 

(-A + ml) wo - g^ {pl + p'J) , 

(-A + m2)a; = g^(j^+j"), 

(-A + m^) po - gp (p^ - p") , 

(-A + m2)p = g^(_;-p_j"), 

-AAo - ep^, 
-AA = ejP, 

where the source scalar and vector densities read 

1=1 

N,Z 



(19) 



p"'^ = E^l'^^ 



and the source currents 



Af,Z 



r'^ = E^»'^^- 



(20) 



(21) 



4=1 



where n andp refer to neutrons and protons, respectively. 
Equations ([TT]) and p^ provide a closed set. Their so- 
lution has to be found iteratively, starting from a rea- 
sonable estimate of the meson fields, the Dirac equation 
dill) is solved yielding the single-particle spinors. From 
the spinors, using (I^IH) and (HU, one obtains the densities 
and currents, which act as sources for the solution of the 
Klein-Gordon equations (fT9|) that provide a new set of 
meson fields. Repeating the procedure until convergence 
results in the self-consistent solution of this set of equa- 
tions (see Ref. [l2] for details). From this solution one 
can calculate physical quantities such as the total energy, 
radii and deformations. 

The actual numerical solution of these coupled set of 
equations is obtained using an oscillator expansion in 



-7.7 
> -7.9 

CD 

E -8.1 F 
H~ -8.3 



? 2.8 

-M 

-: 2.7 

.£ 2.6 



■a^ 




-0.2 h 
-0.4 



-1 1 1 1 1 1 1 1 1 1 1 1 1 — -| 1 1 1 1 r- 



— Self-consistent 
— Frozen Fields 




I I I I I I I I I I I I I I I I I I 1 ^ 



Total 

Neutron 

■ — • — Proton 



I I I I I I I I I I I I I I I I I I I 



Total 

Neutron 

• — ■ — Proton 



_i I I I I I I I I I I I I I I I I I i_ 



0123456789 
B (IQi^ G) 

FIG. 1: (Color online) Binding energy per article, radius and 
/3 deformation dependence on the magnetic field strength for 
^®0. The binding energy per particle also shows the differ- 
ences between using a self-consistent approach and adding 
the magnetic field in a frozen-field configuration on top of the 
bare self-consistent ground-state (see text for details). 



TV = 20 major shells, for which further technical details 
can be found in Ref. [l^l- Details pertaining the new 
terms involved in the inclusion of the coupling to an ex- 
ternal magnetic field can be found in the Appendices A 
and B. 



III. EFFECTS OF THE MAGNETIC FIELD ON 
THE NUCLEAR STRUCTURE 

The effects that the coupling of protons and neutrons 
to an external magnetic field has on the nucleus can be 
classified as: 

Neutron paramagnetism: or Pauli-type magnetism, is 
caused by the interaction of the magnetic field with the 
neutron magnetic dipole moment. It induces a relative 
shift of levels with neutron spins directed along the mag- 
netic field. Since the gyromagnetic factor for neutrons 
is negative (<;„ ~ —3.8263), configurations with the spin 
anti-parallel to the magnetic field are energetically fa- 
vored. 

Proton paramagnetism: as in the case of neutrons, it 
comes from the interaction of the magnetic field with the 
proton magnetic dipole moment. However, the gyromag- 
netic factor for protons is positive {gp = 5.5856), which 
favors configurations where the proton spin is parallel to 
the magnetic field. 

Proton orbital magnetism: or Landau- type magnetic 
response, that couples the orbital motion of protons with 
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FIG. 2: (Color online) Left: neutrons, right: protons. Evo- 
lution of the single-particle levels in ^^0 with increasing ex- 
ternal magnetic field, with frozen nuclear potentials at their 
values for vanishing B. Landau coupling for the protons and 
the coupling of the anomalous magnetic moments is included. 
Blue lines refer to levels with positive angular momentum pro- 
jection Q,, while red lines to levels with negative SI. Solid lines 
indicate positive parity, while dashed lines indicate negative 
parity. The green dots mark the last occupied level, while the 
magenta circle on the proton graph highlights the first level 
crossing at the Fermi energy. 



the magnetic field. It favors configurations where the 
proton angular momentum projection is oriented along 
the direction of the external magnetic field. 

From the single-particle level point of view, there are 
two different effects. The orbital magnetism associated 
with proton ballistic dynamics removes Kramer's degen- 
eracy in angular momentum projection Vi of proton levels, 
and brings those aligned with the magnetic field down in 
energy. On the other hand, the paramagnetic response 
(Pauli magnetism) removes the angular momentum pro- 
jection degeneracy for both protons and neutrons. It is 
thus expected that the magnetic field effect on the single- 
particle structure is more pronounced for protons than 
for neutrons. 

As a first step, it is enlightening to study the effects of 
magnetic fields on the nuclear structure in a simple and 
well known nucleus like ^^O. Fig. [T] shows the evolution 
of the bulk properties of ^^O with increasing magnetic 
field. At first the influence of the external magnetic field 
is counteracted by the currents generated by the breaking 
of time reversal symmetry, including the classical oc B^ 
contribution coming from the orbital coupling. The ra- 
dius of the nucleus and spherical shape show high re- 
silience to the increase in the external magnetic field. 
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FIG. 3: (Color online) Same as Fig. [2l but with a fully self- 
consistent solution of the equations of motion (see text for 
details). 



For field strengths around 5 • 10 ^^ G, there is an abrupt 
decrease in binding energy, associated with increased ra- 
dius and the sudden appearance of oblate deformation for 
the ground state. Such discontinuities are an indication 
that the underlying shell structure has changed in a fun- 
damental way. These jumps in bulk nuclear properties 
can be traced to the single-particle behavior, as can be 
observed in Figs. [2]and[21 they occur when the last occu- 
pied level crosses the first empty level. At that point, for 
even-even nuclei, a reoccupation occurs. A particle is re- 
moved from a level going upwards with increasing B-field 
and brought to a level going downward with increasing 
spin. Since these two levels have opposite angular mo- 
mentum along the symmetry axis, the nucleus becomes 
spin-polarized. Another effect on the structure of nuclei 
which are superfluid for vanishing external magnetic field 
is the gradual disappearance of the neutron and proton 
pairing-gaps with increasing external field. The origi- 
nal shell structure is washed out due to the complicated 
pattern of level crossings, and, as the magnetic field in- 
creases, new magic numbers may appear. 

In the upper panel of Fig. [T]we show two types of cal- 
culations. The full curve correspond to self-consistent 
calculations, where the nuclear potentials change with 
increasing external fields due to polarization of the den- 
sities and due to the polarization currents (nuclear mag- 
netism). The dashed curve (frozen potentials) corre- 
spond to a calculation where, in a first step, the nu- 
clear potentials S'(r) and Vb(r) = 3^^0(1") +3pT'3Po,3(r) + 
eAo (r) are calculated without external magnetic field and 
subsequently these potentials are kept frozen while the 
external magnetic field is switched on. Therefore, in this 



case we have no nuclear magnetism. For small external 
fields, i.e. up to the first level crossing at B fa 5.2 • 10^^ 
G there is practically no difference. Levels with ±51 are 
equally occupied and the corresponding single particle 
wave functions are very similar; their contributions to 
the currents nearly cancel each other and there is prac- 
tically no polarization and no nuclear magnetism. The 
situation changes, however, after the first level crossing. 
Now the nucleus is spin polarized in the self-consistent 
solutions (full curve), we have polarization currents and 
nuclear magnetism, effects neglected in the calculations 
with frozen fields. Thus we observe differences in the 
binding energies and also in the location of the next level 
crossings. 

In order to understand the results of Fig. [1] in more 
detail we consider in the Figs. [5] and [3] the effects on the 
single-particle structure outlined at the beginning of this 
section for frozen fields and for self-consistent fields. Pro- 
ton and neutron level degeneracy is broken in reversed 
directions due to the different sign in their paramagnetic 
behavior interacting with the external magnetic field. 
This degeneracy breaking is more acute in the case of 
protons, where the orbital magnetism plays an impor- 
tant role. Therefore the first level crossing occurs for the 
protons at _B w 5 • 10"^'' G. For frozen fields (Fig. ^ this 
first level crossing for protons has no influence on the 
other proton levels nor on the neutron levels. This is no 
longer the case for the self-consistent solution in Fig. |3l 
where the changes in the nuclear fields caused by polar- 
ization are clearly seen also in the other proton levels and 
due to the proton-neutron interaction also in the neutron 
levels. 

Polarization effects induced by the w and p currents 
due to breaking of time-reversal symmetry are impor- 
tant, and the frozen fields approximation breaks down 
for higher magnetic fields. In the first level crossing one 
proton is removed from the lpi/2 shell and brought to the 
downward sloping orbit of the 1^5/2 shell, i.e. to a level 
with 51 — +5/2. In the second and third panels of Fig. [T] 
we see that the reoccupation corresponds to an increase 
of the proton radius and transition from a spherical shape 
to an oblate deformation. 

Light stable nuclei, like ^^O in Fig. [1] are very stiff in 
their response to the external magnetic field. In this par- 
ticular case the change in binding energy is less than 100 
keV per particle for field strengths less than 5 • 10^^ G. 
The induced currents tend to counteract the effects of the 
magnetic field. The point at which the first level crossing 
occurs may be arbitrarily defined as the minimum field 
strength for which the nuclear structure is significantly 
altered. Thus for ^^O that would be 5 • 10^'^ G, with a 
jump in binding energy per particle of around 400 keV. 
For less stable and/or heavier nuclei, as will be studied 
in the next section, it is expected that this minimum 
field is reduced, due to, mainly, two effects: i) increase in 
the level density around the Fermi energy, and ii) the in- 
crease in the proton orbital magnetic response due to the 
occupation of single-particle orbitals with higher angular 



-8.5 



^^1^i:il->>~4- 






> 

CD 






-9 - 



-9.5 







^ «'Fe 


■ -T- ■ 


5 Fe 


Ti 


..-•-. 


5BNi 


....... «Cr 


.-*-- 


« Zn 



J I L_ 



4 6 

B (1017 G) 



10 



FIG. 4: (Color online) Energy per particle for N=Z nuclei 
around ^ Fe. 



momenta. In particular, this implies that the validity of 
a frozen-field treatment of the coupling to an external 
magnetic is highly dependent on the nucleus under con- 
sideration. In the next section it shall be investigated the 
response of heavier nuclei to the external magnetic field. 



IV. POSSIBLE INFLUENCE ON THE OUTER 
CRUST COMPOSITION 

The outer crust of neutron stars, below the neutron 
drip density, is believed to be composed by well separated 
nuclei positioned on a body centered cubic (bcc) lattice in 
complete thermodynamical equilibrium. Assuming that 
matter in such conditions condenses to a perfect crys- 
tal lattice with a sin gle nuclear species at each site, the 



energy density is 



e = nNWNiA,Z) + e'^ine)+eL{Z,ne) (22) 

where Wn is the mass-energy of the nuclear species, e^ 
is the free energy of the electrons, e^ is the bcc Coulomb 
lattice energy, and tin and rie are the number densities of 
nuclei and electrons. At a given pressure, minimization 
of the Gibbs free energy per nucleon with respect to the 
nuclear species (A,Z) determines the nuclear composition 
of the lattice vertices. 

So far all previous studies (see, for example, [3,27-31]) 
on the influence of strong magnetic fields in the com- 
position of the crust have concentrated on the electron 
part, and assumed that the nuclear binding energy Wn is 
not affected. This is certainly true for weaker magnetic 
fields are below 10^^ G, since such fields do not alter 
the nuclear structure. However, previous studies [8|, |9| 
where the influence of the magnetic field on the shell cor- 
rection energy has been investigated in a simple model, 
found that fields with a strength above that threshold 
might change the nuclear binding energy, and thus mod- 
ify the equilibrium nuclear species on the lattice. In fact, 
changes of a few keV per nucleon might significantly alter 
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FIG. 5: (Color online) Left: neutrons, right: protons. Evo- 
lution of the single-particle levels in *^Fe with the magnetic 
field, including the proton orbital coupling and the anoma- 
lous magnetic moments coupling. Solid (blue) lines refer to 
levels with positive Q, while dashed (red) lines to levels with 
negative Q. Solid lines indicate positive parity, while dashed 
lines indicate negative parity. The (green) dots mark the last 
occupied level. 



the composition, making more neutron-rich nuclei domi- 
nate over the most likely to be found nuclei around the 
^^Fe region ^. 

One important question is the location of this mag- 
netic field strength threshold, in particular for nuclei in 
the vicinity of ^^Fe. Assuming an electron fraction close 
to Ye = 0.5 13^, it can be argued that the most prob- 
able nuclei are those with the same number of protons 
and neutrons. Fig. S] shows the dependence of the nu- 
clear binding energy per nucleon on the magnetic field 
strength for N=Z nuclei for fully self-consistent solutions 
of the RME equations in Eqs. ([TT|) and ([T^ in the re- 
gion close to ^^Fe. For field strengths of 0.5 • lO^^ G 
there are changes of a few tenths of keV in the binding 
energy per nucleon of some species. At higher magnetic 
fields, around 2 • lO^^ G, the hierarchy in binding energy 
of the most bound nuclei is altered and thus it is expected 
that the composition of the outer crust is substantially 
altered. 

As in the example with ^^O, the origin of these discon- 
tinuities in the binding energy per particle can be traced 
back to the single-particle structure. Fig. [5] shows the 
single-particle level spectra with respect to the external 
field for ^®Fe. The diagram is very similar to that of ^^O, 
only the larger level density around the Fermi energy 
and smaller single particle gap reduce the minimum field 
for which the magnetic field produces structural changes. 



Similar diagrams are found for all of the other nuclei pre- 
sented in Fig. [5l and in those in the vicinity of ^^Fe. Or 
course, details pertaining the minimum field that induces 
a different single-particle level occupation scheme depend 
very much on each particular nucleus. 

Concerning the possible changes in the hierarchy of 
binding energy per particle with increasing magnetic 
field, it has been found in this work that several nuclei 
(e.g. ^^Fe and ^^Fe) overbind ^^Fe for extended ranges 
of external magnetic field strengths. In that regard, the 
frozen field solutions show a different behavior than the 
self-consistent ones. The field ranges for which this hier- 
archy changes occur are different and the magnitude is 
typically off by a couple hundred keV as compared with 
the fully self-consistent solutions. This ordering accord- 
ing to binding energy is one of the factors that influences 
the final composition of the outer crust in neutron stars, 
and quantitative predictions should be done using the full 
self-consistent formulation. 

It is also interesting to study the minimum field (de- 
fined as the field at which the first level crossing occurs) 
for isotopic chains in the iron region, since it provides 
an indication of the possible effects on the neutron star 
composition. Fig. |6] shows the magnetic field at which 
the first level crossing occurs, for different isotopic chains 
close to iron. It provides an intuitive idea of how intense 
the magnetic fields have to be on average in order to af- 
fect the nuclear structure significantly. For isotopes close 
to ^^Fe this value is approximately between 0.5 • 10^^ G 
and 3 • 10^^ G. For heavier as well as neutron-rich nu- 
clei, which exist at higher densities in the crust, a sharp 
decrease of this minimum field is expected as mentioned 
previously. Thus, for fields around 10^^ G it is reason- 
able to expect changes in the crust composition. How- 
ever, and because of the strong dependence of this min- 
imum field on the nuclear species, it is not possible to 
predict the effect on the composition without actually 
performing the minimization of Eq. (j22p . This calcula- 
tion is, however, computationally very demanding and 
would also require refinements in the model like the in- 
clusion of pairing (even though it is reasonable to expect 
that magnetic fields damp and eventually cancel it) or 
the proper inclusion of one-pion exchange terms in the 
effective Lagrangian. Therefore, it is well outside of the 
scope of the present work. 



V. CONCLUSIONS 

The infiuence of strong magnetic fields on nuclear 
structure has been studied using a fully self-consistent 
covariant density functional. It has been found that a 
field strength of at least 10^^ G is needed to appreciably 
modify the nuclear ground state. For sufficiently high 
fields these effects cannot be studied using a frozen-field 
approach since there are level rearrangements, causing 
spin-polarization and induced currents and thus a self- 
consistent model is required. It is the advantage of co- 
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FIG. 6: (Color online) Minimum magnetic field for which 
the first level crossing at the Fermi energy occurs for isotopic 
chains around ^^Fe. 



variant theories that the these currents can be taken into 
account without any additional parameters. The mini- 
mum magnetic field that changes the nature of the nu- 
clear ground-state is very much dependent on the nu- 
cleus, and its effects on the nuclear binding energy per 
nucleon range from a few tenths of keV for B sa 0.5T0^^G 
to a few hundreds of keV for the maximum field theoret- 
ically possible B w 10^^ G. 

No neutron star has yet been observed with such in- 
tense magnetic fields B > 10^^ G, even though theoreti- 
cal models hint that such objects exist [lOjllll- In such 
a case, the composition of its outer crust might be rad- 
ically different from that of normal neutron stars. The 
relevance of this change in composition depends on the 
abundance of normal neutron stars compared with that 
of magnetars. Until now, observations suggest that mag- 
netars are not so common in the universe and thus it is 
unlikely that the magnetic field effects on nuclear struc- 
ture play an important role in global astrophysical ob- 
servables like element abundances. However, changes in 
the composition of magnetars might be relevant in the 
study of different phenomena in these particular kind of 
neutron stars, for example elastic properties of the crust 
3J], pulsar glitches [35|, or cooling [36]. With the inclu- 
sion of a proper pairing interaction, the covariant DFT 
model presented in this work can be used to perform a 
quantitative exploration of all these questions. A sys- 
tematic study of changes in the composition of the outer 
crust will be presented in a upcoming publication. 
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VI. APPENDIX A: OSCILLATOR MATRIX 

ELEMENTS FOR THE PROTON ORBITAL 

COUPLING 

A. Dirac equation 

The Dirac equation (jlip together with the meson field 
equations p^ are easily solved in a harmonic oscillator 
basis. The procedure is described in great detail in Ref. 
[12| . The only new term coming in the solution of the 
Dirac equation is a • V where V is given in Eq. ([T5|). 
Using the set of a matrices in the spherical tensor basis, 
(a+, a^, as) with 



cr+ \ f (7- 

"+ = i a+ ) ' "- = U- 



allows to write 

aV = a+y- + a^V+ + asV^, 
with 



(23) 



(24) 



(25) 



and the internal self-consistent magnetic potential V 
given by the solution of the Klein-Gordon equations. Us- 
ing an oscillator expansion for the spinors 



'ipi = 



Enfn\n) 



(26) 



where \n) are the axially symmetric harmonic oscillator 
wave functions, determined by the quantum numbers n = 
{nz,nr,mi,ms). 



\n) = \n^nrmims) = 0„Jz)0|,"'l(r)^==Xm.('S) (27) 

V 27r 



with VL — mi-\- TTig . The term ol-V in the Dirac equation 
is then in matrix form 



Bnn' \ f In 
Bn'n y I ign' 



(28) 



with 



Bnn' = -~{n\V-c7+\n') - {n\V+cj^\n'). (29) 

The oscillator matrix elements can be written as 

{n\V~(J+\n') ^+6rn[,m^+i{n^n^mi\V^\n'^ji'^m'i) (30) 
(n|F+a_|n') =-<5™;,™,-i(n,n,m,|F+|«m;). (31) 

B. Currents in coordinate space 

To solve the Klein-Gordon equations, the source terms 
have to be transformed from oscillator to coordinate 



space. Expressions for the scalar and vector densities 
are given in Ref. [12l |. The currents are defined as 



j(^) ^^i^\{r)oLikt{r) 



(32) 



Because of axial symmetry the third component vanishes, 
and we have to consider only two of the components: 

j+{r) = XI ^1^+^^ ^-{r) = X! V-la-V-j- (33) 



Since a+ = (q;_)^ we have j_ = j^, so only one of them 
needs to be calculated explicitly. With the Dirac spinors 
in Eq. ^ we find 






(34) 
(35) 



with 



j{r, z)=J2Yl fn9n'^n{r, z)$„'(r, z)5„,,^„,^+, (36) 

i nn' 
i nn' 

where the $„(r, z) are the oscillator wave functions with- 
out spin dependence, i.e. ^ni^jz) — </>«, (2;)0n^ (r). 



C. Klein-Gordon equation oscillator matrix 
elements for the vector terms 



In the spherical tensor basis, the K-G equations read 

{^A + m^)w' =gf ^ = +,-,3 (37) 

for the massive mesons. The functional form of the cur- 
rents, j±{r,9,z) = ±ij(r, z)e^*^, suggest the following 
ansatz for the potentials: w±{r,9,z) = ±i w(r, z)e^'^. 
Inserting both in the K-G eqs, and eliminating the angu- 
lar dependence in 9 



dl - -dr -dl + ^+fi?) w(r, z) = j(r, z) (38) 



The oscillator matrix elements for the Laplacian can be 
found in [l^]. It is only a matter of including the oscil- 
lator matrix elements for the 1/r^ term, which can be 
accomplished trivially in coordinate space. 



VII. APPENDIX B: OSCILLATOR MATRIX 

ELEMENTS FOR THE COUPLING WITH 

INTRINSIC MAGNETIC MOMENTS 

The coupling of the magnetic field with the anomalous 
magnetic moments of protons and neutrons introduces a 
new term in the single-particle Dirac equation dill) with 



AT3 — '^TaH'N r^fJfil^^ 



(39) 



For a constant magnetic field B of the form ((T6| it is 
easy to show that 



2^" 



(T 

a- 



B= S B, 



and for B — (0.0, B) we have 



'^"^ - -KrsflN'^sB, 



(40) 



(41) 



Finally, the oscillator matrix elements can be written as 



{n\ixil'\n') =±S„^n'Jn^n'Jmim[Sm,Tn'^Kr3lJ-NB (42) 



for nis = ±i. 
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